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Abstract 

In this paper we discuss the ergodic properties of quantum conservative systems by 
analyzing the behavior of two different models. Despite their intrinsic differencies they 
both show localization effects in analogy to the dynamical localization found in Kicked 
Rotator. 



1 Introduction 

The study of quantum systems which are classically chaotic is known under the 
general term of Quantum Chaos. This subject is very interesting and quite new. 
Indeed, while the statistical properties of chaotic classical systems are well described 
by the ergodic theory 0, much less is known about Quantum Chaos. 

Since the preliminary investigations in the 70 's |2| it became quite clear that 
one of the main characteristics of the classical chaotic motion, namely the diffusion 
in the phase space, is suppressed by quantum mechanics. Despite the fact that this 
phenomenon was discovered in a particular time-periodic Hamiltonian system, (the 
so-called Kicked Rotator model) it was later predicted and observed in other, more 
physical systems, like the Hydrogen atom in a microwave field|3|]. The quantum 
suppression of classical excitation has been called "dynamical localization" || . 



On the other hand, bounded conservative systems, were analyzed from the point 
of view of the statistical properties of eigenvalues and eigenvectors Important 
analogies were established between the statistical properties of Hamiltonian spec- 
tra and those derived from random matrices for which a fairly deep mathematical 
analysis had been developed^]. In particular the spacings distribution for near- 
est neighboring levels of a classically chaotic system was found to be close to that 
derived from random matrices belonging to a given symmetry class || (Gaussian 
Orthogonal (GOE) or Gaussian Unitary Ensembles (GUE) depending on whether 
the Hamiltonian is invariant or not under time reversal). Other approaches are 
based on nonlinear a-models and supersymmetric techniques. 

In this paper we investigate two simple conservative systems in order to show 
that the dynamical localization can appear in this case also. Roughly speaking, in 
an ergodic conservative system the motion takes place on the whole energy surface. 
Correspondingly, in the quantum case, an eigenfunction can be extended or localized 
inside the energy shell. 

A main difficulty in understanding quantum chaotic motion arises from the fact 
that a bounded conservative quantum system is characterized by a discrete energy 
spectrum with a quasi-periodic motion and this can be hardly compared with the 
typical features of the corresponding classical chaotic motion which is character- 
ized by a continuous spectrum. At first glance this would seem a failure of the 
correspondence principle. However the distinction between discrete and continuous 
spectrum becomes meaningful only on infinite time. On finite time scales, instead, 
the quantum motion can be chaotic as the classical one. 

A discussion of time scales, which play a fundamental role in the study of Quan- 
tum Chaos, is given in Section 2. In Section 3 we will analyze a model of classically 
ergodic conservative system, the so-called Wigner Banded Random Matrix ensem- 
ble, where the dynamical localization phenomenon is shown to take place. Finally, 
in Section 4, we will investigate the problem of localization in a more physical 
dynamical system : the Bunimovich stadium. 

2 Time Scales 

What are the physically meaningful time scales in quantum dynamics? To answer 
the question let us compare the quantum and classical evolution starting from the 
same initial conditions. Namely, let us consider the evolution of a quantum narrow 
packet and of a beam of classical orbits initially centered in the same small area 
of phase space. As it is known, the narrow packet will follow the beam of classical 
orbits as long as the beam remains narrow. Due to the exponential instability of 
classical chaotic motion, this can happen only up to a time tg, called Ehrenfest 
time, 

t E oc log(I/h) 
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where / is a typical value of the action of the system. This time scale, first intro- 
duced in J7|, is very short. Yet, according to the correspondence principle, it grows 
to infinite as h — > 0. For times t < tE the quantum dynamics is chaotic as the 
classical one, including the exponential instability of the motion. 

Another, more interesting and longer time scale, is related to the Heisenberg 
uncertainty principle. Indeed, if p* is the level density of the operative eigenstates 
involved in the quantum dynamics (namely the eigenstates which enter in the ini- 
tial condition) then the discreteness of the spectrum can be resolved by quantum 
dynamics only after a time £b[Q : 

t B ^ hp* (1) 

This time scale is related to the process of relaxation to the quantum steady 
state. Indeed, for times t < £b, quantum and classical averages of dynamical 
observables will be close to each other while, for t > ie , the quantum system will 
"see" the discrete nature of the spectrum and consequently it will reach a stationary 
regime. Quite obviously, as required by the correspondence principle, this time also 
goes to infinite when h — > 0. Notice, however, that as h — > 0, the time tB diverges 
to infinite according to a power law dependence. Therefore the quantum relaxation 
time is much larger than the time tE which characterizes the quantum instability 
of the motion. This means that the quantum diffusion and the relaxation process 
take place in absence of exponential instability, which is confined in the small time 
interval tg. 

The nature of the quantum steady state will depend on the comparison between 
the time tB and the ergodic time t erg for the classical relaxation to equilibrium. 
Indeed, if L is a typical length scale over which classical diffusive motion takes 
place with diffusion rate D, then the ergodic time t erg can be estimated as: 

t erg ~ L 2 /D (2) 

It is clear that, if the quantum relaxation time is is larger than the classical 
one t erg , then the quantum steady state will be close (apart quantum fluctuations) 
to the classical one, given by the microcanonical ensemble. On the other hand, a 
more interesting situation appears when ts < t erg . In such a case, the quantum 
distribution will relax to a steady state which is not ergodic but localized around 
the initial excitation. 

In conclusion, the nature of the quantum steady state is determined by the 
parameter A, which has been called "ergodicity parameter" : 

A 2 = tB/t erg (3) 

Systems characterized by A ^> 1 will relax toward a quantum ergodic steady 
state, while, for A < 1, the quantum stationary state is localized. In the follow- 
ing sections we will discuss, on specific examples, the mechanism through which 
quantum dynamical localization can actually take place in conservative systems. 
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3 The Wigner Band Random Matrices Ensemble 



Apart from few recent papers |§, the main interest in the Random Matrix Theory 
has been related with full random matrices ||. However, for real quantum Hamil- 
tonian systems, matrix elements decay on moving away from the main diagonal 
and this lead to the consideration of matrices with a band structure, in which ma- 
trix elements are different from zero only inside a band of size b around the main 
diagonal. 

Let us consider, therefore, the following band random matrix(WBRM) : 

H m ,n ^n^m,n ^m,n (4) 

which can be taken to describe an Hamiltonian system of the type H = Hq + V 
where Hq is integrable and V is a perturbation which renders the Hamiltonian H 
ergodic and mixing. The unperturbed energies e n are assumed to be distributed 
according to a Poisson law. The off-diagonal matrix elements v m ^ n are taken as 
Gaussian random numbers with zero average and variance v inside a band of size b: 

{v^ nn ) = v 2 for \m — n\ < b. (5) 

Outside the band b, the matrix elements are zero. This model is then defined 
by three parameters : b, v and the average density p of unperturbed levels: 

p = 1/ (e n - e n _i). (6) 

This model, introduced by Wigner M, has nowadays attracted much attention 



[[10| , |LT| , and several interesting results have been obtained. Here, we only briefly 
mention those connected to the localization problem. 

The quantities of interest are the Strength Function or Local Density of States 
(LDOS): 

w(E\E ) = Y.(\ME m )\ 2 )n 6{E - E m ) (7) 

m 

and the averaged eigenfunctions distribution in the energy space: 
W(E \E) = J2(\ME m )\ 2 ) m S(E - El) (8) 

n 

Here, ip n (E m ) is the n-th component of the eigenfunction having E m as eigen- 
value and the averages ((. . .)„, (. . .) m ) have been done respectively over a small 
number of n values close to Eq and over those eigenfunctions having an eigenvalue 
E m close to E. 
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The non-perturbative regime is defined by the condition that the strength of 
the perturbation must be larger than the average unperturbed levels spacing : 



pv > 1. 

Above the perturbative regime, the shape of LDOS depends on the Wigner 
parameter: 

pv 2 

q = Wp 

Here the numerator represents the spreading width induced by the perturbation 
(given by the Fermi golden rule), while the denominator stands for the width, in 
energy, of the band matrix. 

It is possible to show that, when g > 1, the LDOS has a semicircle form with 
diameter AE = AvV2b (semicircle regime (SC) ). In the other regime, characterized 
by q <C 1, the LDOS has a Lorentzian shape with the main part inside a width 
r = 2npv 2 « AE = 2b/ p (Breit-Wigner regime (BW) ). 

Due to energy conservation, the number of states occupied by an eigenfunction 
will be bounded from above by the ergodic localization length l e ~ pAE which 
gives the maximum number of states which can be coupled by the perturbation. 
We then consider, as localized, those states with a localization length significantly 
less than the maximum one £ e , namely we consider localization inside the energy 
shell AE. As a measure of the degree of localization we can take, for instance, the 
inverse participation ratio, 

£=l/Y,\ME m )\* (9) 

n 

It can be shown [fLT|l that the ratio j3i oc = £/£ e obeys the following scaling relation: 

Pioc = y ~ 1 - e" A (10) 

where the parameter A, which plays the role here of the "ergodicity parameter", 
is defined as 

I b 2 



In Eq. (|lT|) , £oo is the localization length for band random matrices with infinite 
unperturbed density p(T^, [13[. The parameter A is called ergodicity parameter since, 
when A > 1 , one has (3 ~ 1 and therefore i ~ t e . On the other side when 
< A < 1 one has (5 ~ A < 1 and t < t e . 

The global properties of the eigenfunctions are also connected with the statistical 
properties of the spectrum. Indeed, it was found that effects of localization manifest 
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in the repulsion of neighboring levels. To be more precise, one can show||14j| that, in 



case of localization, the distribution of neighboring levels spacing obeys the following 
relation: 

p(s) = As p exp [ - n 2 f3s 2 /lQ - (B - n/3/4) s] (12) 

where A, B are constants obtained from normalization. An interesting fact, which 
so far has not yet received a theoretical explanation, is that the numerical value 
of the repulsion parameter (3 in (12) turns out to be very close to the localization 



parameter Pi oc [ll 



4 The Bunimovich Stadium 

In the previous section we have shown that quantum dynamical localization can take 
place in the model of WBRM. Even though WBRM are believed to describe the 
qualitative properties of conservative, classically chaotic, Hamiltonian systems, it is 
highly desirable to analyze a more realistic Hamiltonian system. To this end, two- 
dimensional billiards are very convenient objects to study since they have very clean 
mathematical properties, from complete integrability (e.g. the circle) to complete 
chaotic motion (e.g a dispersive billiard). Moreover their classical and quantum 
dynamics can be numerically studied with sufficient good accuracy. Finally, modern 
laboratory techniques allow for quite accurate experimental investigations. 

In the following we consider a well known billiard model, the so-called Buni- 
movich Stadium, which consists of two semicircles, with radius R — 1, connected 
by two straight lines with length 2a. Inside this bounded two-dimensional region 
we consider the motion of a point particle with mass m — 1 colliding elastically 
with the boundary. The classical dynamics depends only on the ratio e = a/R and 
it can be rigorously proven to be ergodic and mixing for any e ^ 0. When e ~ 1 the 
relaxation time to statistical equilibrium is very short, just few collisions with the 
boundary. Here we are interested to the case e < 1 when the stadium is very close 
to the circle. For the billiard in a circle, there exist two constants of motion, the 
energy E = mv 2 /2 and the angular momentum I = l z k (here k is the unit vector 
perpendicular to the plane of the stadium and identifying the z axis). For e > the 
rotational symmetry around the z-axis is broken and l z is not a constant of motion 
any more. Nevertheless, if e is sufficiently small, one can expect that the angular 
momentum will change slowly in time. Indeed, the angular momentum can vary in 
the interval \l z \ < l max = m(R + a)\v\ ~ \ / 2mR 2 E; now, if for example at t — 
the angular momentum is zero, then it will evolve in time in a diffusive way, with 
diffusion coefficient D, until the system will reach the equilibrium state given by 
the microcanonical ensemble. The ergodic time t erg , namely the time necessary to 
reach the statistical equilibrium, can be estimated as: 

terg * tjD (13) 
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The diffusion coefficient D can be computed numerically ||15|| . In terms of the 
rescaled angular momentum L = l z /l ma x, the diffusion coefficient, in the number 
of collisions with the boundary, is given by D = 1.5 e 5 / 2 . By taking into account 
that the time interval between two successive collisions is t c ~ E^ 1 ^ 2 one has that, 
neglecting numerical constants (m = R = 1) , 

ter 9 ~ ^t e ~ e^E^ 2 (14) 

The classical dynamics of the billiard can be approximated by the following 
area-preserving map which gives the change of L and of its related conjugated 
variable 9, between two successive collisions with the boundary ||15||: 

L = L + e sin 6 sgn(cos 6) 

= # + 7r-2arcsinZ (15) 



Here e = — 2 e sgn(L ) y 1 — L\ and L is the initial value of the angular momen- 
tum. 

The map fllS]) represents a first order approximation to the real dynamics when 
the initial rescaled angular momentum L is not too large. It also give rise, for 
s « 1, to a diffusive motion with a diffusion rate D oc e 5 ^ 2 in agreement with 
numerical computations on the real billiard. 

The map description allows to understand the properties of the quantum dy- 
namics in analogy to the Kicked Rotator model0. Indeed we know that, above 
the quantum perturbative regime e > H, quantum dynamics will follow the classical 
diffusive motion up to a certain time, called break-time, (see also Eq.(0)) given 
by 75 ~ -Do 05 where r& is measured in the number of collisions and Dq is the 
dimensionless diffusion coefficient. In physical units, this time is given by: 

t B ~ e^E^/h 2 (16) 
This allows to estimate the ergodicity parameter A (Eq.@) 

A 2 = t B /t erg = Ee 5 /h 2 (17) 
The value A = 1 gives the critical energy value: 

E erg ~ U 2 /e 5 (18) 

above which we expect quantum ergodic behavior. On the contrary, for energies 
E < E erg , we expect dynamical localization. Since the total average number of 
states up to the energy D0 is given by 

(M(E)) « ^ 2 E (19) 



7 



where A is the area of the billiard, one can estimate the level number N erg ~ e -5 
above which quantum ergodic behavior is expected. Numerical computations show 
that, in agreement with the above estimates, only above N erg the levels spacing 
distribution is very close to a Wigner-Dyson distribution, characteristic of the GOE 
ensemble. 

One can also numerically compute the repulsion parameter /3, using Eq . (|T2|) and 



one finds 16 



P ~ 1 - e" A (20) 



where A is the ergodicity parameter ([TTp. This quite surprising result reinforces the 
similarity between WBRM and the quantum behavior of classically chaotic systems. 

Finally , we would like to notice that the ergodicity parameter A turns out be 
proportional to the dimensionless conductance: 

PD 

where E c is the Thouless energy and A is the average levels spacing. Indeed 
taking into account that 

E c ~ -A (22) 

^erg 

and 

A ~ A (23) 
where — h^^- ~ h^ 1 (see Eq.(|T0)) is the Heisenberg time; one gets 



Q = IE. 111 oc A (24) 

5 Conclusions 

In this paper we have discussed the problem of localization in conservative, classi- 
cally chaotic, Hamiltonian systems. The existence of localization in such systems 
would restrict quantum distributions to smaller regions of phase space than classi- 
cally allowed, and would therefore introduce significant deviations from ergodicity. 
As surmised in ||17|1 , this lack of quantum ergodicity may lead to interesting conse- 



quences for quantum equilibrium statistical distributions. 
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